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ABSTRACT 

We consider the post-Newtonian approximation for the Dedekind ehipsoids in the case of ax- 
isymmetry. The approach taken by Chandrasekhar & Elbert (1974, 1978) excludes the possibility 
of finding a uniformly rotating (deformed) spheroid in the axially symmetric limit, though the 
solution exists at the point of axisymmetry. We consider an extension to their work that permits 
the possibility of such a limit. 

Subject headings: Dedekind ellipsoids, non-axisymmetric, stationary 



1. Introduction 

The Dedekind tri-axial ellipsoids are an ex- 
ample of non-axisymmetric, but stationary solu- 
tions within Newtonian gravity. Due to internal 
motions, they are, in fact, stationary in an in- 
ertial frame. When addressing the question of 
whether or not stationary, but non-axisymmetric 
solutions are possible within General Relativity, 
this property makes the Dedekind ellipsoids a nat- 
ural choice upon which to base one's considera- 
tions. It was, in part, with this question in mind 
that Chandrasekhar & Elbert (1974, 1978) turned 
their attentions to the post-Newtonian (PN) ap- 
proximation of the Dedekind ellipsoids. In a pa- 
per from the same series, Chandrasekhar (1967b) 
had already considered the axisymmetric limit of 
the PN Jacobi ellipsoids at length and was able to 
show that it coincides with a certain PN Maclau- 
rin spheroid (just as their Newtonian counterparts 
coincide at the point of bifurcation). This is re- 
lated to the fact that the PN figures were chosen 
to rotate uniformly. On the other hand, the PN 
velocity field chosen in Chandrasekhar & Elbert 



(1978) excludes the possibility of uniform rotation 
in the axisymmetric limit although it is possible in 
the axisymmetric case. This restriction seems nei- 
ther natural nor advisable in the context of trying 
to settle the question as to the existence of rel- 
ativistic, non-axisymmetric, stationary solutions. 
The nai've expectation is that the axisymmetric 
PN Dedekind ellipsoids contain the PN Maclaurin 
spheroids in the axisymmetric limit (up to arbi- 
trary order). 

In this article, we begin in § 2 by examining the 
axisymmetric case of a generalization to the solu- 
tion presented in Chandrasekhar & Elbert (1978). 
We proceed in § 3 to consider a (continuous) limit 
to axisymmetry. In § 4 the connection to the post- 
Newtonian Maclaurin spheroids is examined. 

2. The Axisymmetric Solution of a Gen- 
eralization to Chandrasekhar and El- 
bert's Paper 

We consider a generalization of the PN Dedekind 
ellipsoids presented in Chandrasekhar & Elbert 
(1978) (referred to from here on in as CE78) in 
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which wc add post-Newtonian terms to the veloc- 
ity. We comply with the notation used in CE78 
and refer the reader to the definitions there for 
the various quantities. The post-Newtonian con- 
tributions to the velocity, which we introduce here 
are 

Svi = alwiX2 + {qi + q)x\x2 + rix\ + tiX2xl 
Sv2 = alw2Xi + {q2 — q)x\x% + r2x\ + t2X-i_x\ 

SV3 = q3XiX2X3, 

(1) 

where the terms with wi and W2 have been added 
for reasons that will be made clear when we dis- 
cuss the solution. Note that we could eliminate 
one constant by introducing variables to denote 
qi + q and 92 — q, but choose instead to retain the 
notation in CE78.^ 

The Newtonian ellipsoid is characterized by the 
semi-major axes oi > 02 > 03. Let us assume for 
the moment that, as in the Newtonian setting, the 
axisymmetric case is obtained by considering 02 = 
ai , an assumption that will be verified shortly. In 
this case, the index '2' in the index symbols Aijk... 
and Bijk... used in CE78 and discussed at length 
in § 21 of Chandrasekhar (1987) can be replaced 
by '1' as is evident from the definitions. Using 
the relations given in that book, it is possible to 
reduce all the index symbols to Ax and A2. At 
the point 02 = ai, the value for Ai (and thus A2) 
is given by (36) in § 17 of Chandrasekhar (1987). 
Furthermore, (2) from Chandrasekhar & Elbert 
(1974) shows us that 



Q2 = —Qi, 



(2) 



where we define the symbol = to mean that the 
expression is evaluated at the point 02 = ai , i.e. 

CL,=„, = ^L,=„, is denoted hy C ^ D. (3) 
The value for as can be found from the equation 

alalAi2 = alAs, (4) 

which holds for the Dedekind (and Jacobi) ellip- 
soids, and gives the value 



— =0.5827241661.... 

ai 



(5) 



-'The three- velocity t)' in CE78 does not refer to the spa- 
tial components of the four-velocity u°' = dx°'/dr, but is 
instead defined as = dx'/dt = u'/m". 



Throughout this paper, a^, is to be understood as 
a function of ai and 02, given by (4). 

We can now consider the integrability condi- 
tions for the pressure and the continuity equation. 
We again follow CE78 and shall refer to the equa- 
tion numbers there by adding a prime. It turns 
out that (38') (of CE78) remains unchanged de- 
spite the modification to the velocity, so that we 
find 



a ^ 

93 = 

and then from (24') that 

a 

92 = -91- 



(6) 



(7) 



Equation (28') is identically fulfilled for 02 = ai, 
moaning that qi is left undetermined, in contrast 

to the general case. 

With the changes to the velocity, equations 
(30') and (31') gain the additional terms {aiQ2Wi + 
alQiW2)xi and {alQ2'Wi + alQiW2)x2 respec- 
tively. Equations (32')-(38') remain unchanged. 
Equation (32') yields 



r2 



-n 



and (37') gives 



t2=h 



(8) 



(9) 



(we shall see shortly that each ti becomes zero). 
There are additional terms in (39') corresponding 
to adding —{aiQ2Wi + a\Q\W2)l'2. — o|(5i(z/;i — 
«;2)/2 both to^ and olT^ . 

Requiring for the new velocity that its normal 
component vanish on the surface leads to a change 
in (50') and thus the resulting equations (52')- 
(56') by which the terms with S1—S2 are modified. 
They now become 



Si — S2 — 



Wi + W2 



2Qi 



(10) 



Using equations (2), (7) and (8), we can subtract 
equation (54') from (55') in CE78 to arrive at 

q + qi-n = lQi{4S3 + S^). (11) 



^We use the superscripts '67' and '78' to distinguish the 
quantities defined in Chandrasekhar (1967b) from those in 
Chandrasekhar & Elbert (1978). 
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Next wc turn our attention to the system of 
equations (58') from CE78^. In the case being 
considered here, the first of these equations be- 
comes 

= qJ^ h 4Q22 2 '3i2 



2al 
3^' 



2al 



^^3(453 + 54) 



i=3 



2ai 



= |Q?(453 + Si) - ^^3(453 + 54) 
+ ia?^iiii(4S3 + 54), 



(12) 



where the values for the a's and their axisymmet- 
ric Umits can be found in Appendix A, the u's 
are given in Chandrasekhar (1967b) (C67b)^ equa- 
tions (72) and (73) and where we made use of (11) 
from the current paper. The unique solution to 
this equation is 



54 = -453 



(13) 



as it is for the analogue equation (100) of C67b 
despite the fact that the term with q~~ri is absent 
there. With the result (13), (11) becomes 



The third minus the second of equations (58') 
is the analogue of equation (101) in C67b and is 
in fact precisely the same equation despite the dif- 
ferent definitions for Pi 



2 



— Qu " ~4Q22 ~2Q 



111 



^2"-3 /-i78 

4 ^;23 



± p78 p78 "3 I d78 D78^ 
— -Til — fo- 



22 



"dfpTs p78\ 
~ TIV-'^Sl ~ ■'^23 ) 



..78 



= an-4|--|(4f-4|) 



-24As{'ialSs + alS,) 



(i) 
^22 



(i) 
'31 



= -2^A3(^a?53 + a255) 

+ (a? Am + 6a^Aii3 - 70^^1113) 
X {^a\S^ + a\S^) . 

The unique solution to this equation is 



S5 



3al 



■S3. 



A3 



(16) 



(17) 



We can use (56') together with (9), (13), (15) 
and (17) to conclude that 



q + qi-ri=0 



(14) 



and for equation (53') from CE78^, or equivalently 
the sum of (54') and (55'), we find 



Dl — D2 = ~3'-'3- 



(15) 



^Please note that wc have been unable to reproduce the 
values from Tabic 1' in CE78 that result from solving (58'). 
A detailed discussion can be found in Appendix A. 

*As mentioned in C67b, the u's belonging to the dis- 
placements ^^^^ and are generated by cyclically per- 
muting the indices. The precise meaning is best under- 
stood via the example that n^*^ can be generated from 



„(3) 



ia|(a|_B23l - aiB22l)- 



aj(a2-Bi23 — aJSiis) and becomes 



(4) 



® In (53') of CE78, the factor Qi is missing from the term 
with (Si - 52). 



t2 = 0. 



(18) 



Equation (47') of CE78 tells us that the bound- 
ing surface is axisymmetric to the first PN order 
if and only if (13), (15) and (17) hold. The PN 
velocity field of CE78 can then be seen to be ax- 
isymmetric in the limit we are discussing, when we 
additionally require 



W2 



-Wi. 



(19) 



Using what has been shown above, the third 
equation of (58') in CE78 can be used to find the 

value of 5*3 (where the relationship between the 
a''^'s and the a^^'s can be found in Appendix A.l) 



= alQll 



4^78 
-'3V33 



2„2^78 



ai a' 



1 "^3*^31 



130 



ii _4 p78 I „4 p78 „2^2 p7« I -™ _2„2 J c 
- (^1^11 + «3-n33 ~ «llC'3-^31 + -^^103^353 
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= atall + atall-alaW3l 

a^Qi / N 130 9 o , „ 
+ ^iq + ri) + —alalAsS3 



i=3 

- „4 67 , 4 67 _ „2 2 67 



2 2 



+ 



+ 



afQi 



130 2 2 . ^ 
2 ''I + -^aifla^s^'a 



(20) 







i=3 



The solution for S3 is then 

53 !^ -0.01742648312 + 0.1061462885 n, (21) 

the analytic expression of which can be found in 
Appendix (B). 

We now turn to the fifth of equations (58') to 
solve for Si . The equation is 

= atQll-atQll + alQl'-alQl' 



^s-" 33 



+ 2alal {3Si + fS^ 



4 67 4 fi7 , 2 fi7 2 fl7 ^fQl 

= 441 - ai«u + 44 - 44 - 

+ 2al4 (3^1 + f 53) - 4QiWi 



i=l 
5 



+ 



■ ^ Si {a\uil - 4'4il + 4'4i - 4'4i^ 

1=3 

(22) 

and the solution is 

Si w -(0.2836731908 + 0.7419729757 ri 



+ 1.121542227 wi), 



(23) 



cf. Appendix (B) for the analytic expression. The 
fourth equation is then identically fulfilled. We 
have obtained a solution to all the equations at 
the point 02 = fli and have two remaining con- 
stants, wi and ri (although q and qi are not de- 
termined, they always appear in the combination 
q + qi, which is equal to ri, cf. (14)). 



3. The Axisymmetric Limit of a General- 
ization to Chandrasekhair and Elbert's 
Paper 

Before discussing the solution obtained above, 
we consider the solution to the post-Newtonian 
equations not at the point 02 = ai, but in the 
limit 02 — > ai. The equations listed above, are 
also obtained as limiting relations. However, in 
the limit, we also obtain two new equations, one 
of which allows us to determine lima^^ai Qi and 
the other, say lima2_^a^ ri. 

Equations (24'), (28') and (38') of CE78 provide 
a system of three linear equations for the quanti- 
ties qi, q2 and 93. After solving this linear system, 
the limit 02 — >■ ai can be taken to give 

qi -6^/2B^l (^^4Biii + J^ns) 



-(2e2 + l)'03 



2.827158725, 



(24) 



where we have defined the eccentricity 



e := 



\/^~4/4 



(25) 



and where the explicit expression for Qi is 



lim Qi= lim _ / 86^(1 (26) 

02^-ai a2^-ai V 3 + 86^ - 8e^ ^ ' 

(we remind the reader that 03 depends on ai and 

a2 via (4)). 

The foiuth of equations (58') is identically ful- 
filled for a2 — ai. Therefore, we introduce 



e := 



2 / 2 



(27) 



and expand the quantities involved and solve to 
first order in s to give 

lim n w 0.02880590648- 1.75 lim wi. (28) 

02— ^ai 02— ^ai 

We provide the analytic expression in Appendix B. 

Strictly speaking, we have to show that the 
fourth of equations (58') is fulfilled to all orders 
in s to be certain that (28) is continuously con- 
nected to the PN Dedekind solutions. We were 
able to solve the whole system of equations along 
the post-Newtonian Dedekind sequence for arbi- 
trary wi and W2, meaning that the limit presented 
here can be tacked on continuously. 
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4. Discussion 

The axially symmetric post-Newtonian solu- 
tions we have generated depend on two parameters 
or one if we require tliat the solution be continu- 
ously connected to the PN Dedekind 'ellipsoids' 
with the velocity field (1). The solutions are not 
uniformly rotating in general. If we add this con- 
straint, then requiring that the four- velocity be 
shear-free tells us that 



ri — >^ (shear free) 



(29) 



must hold. 

We now show that with this additional con- 
straint, the solution is indeed the PN Maclau- 
rin solution (thereby demonstrating that the shear 
free condition is not only necessary, but also suf- 
ficient for uniform rotation in our case). Let us 
first note that upon taking into account the results 
above and in particular Q2 — — Qi, the compo- 
nents of the velocity become 



W3 -0. 



Wi X2 



Wi Xi 



(30) 



This is precisely the form of the velocity for the 
post-Newtonian Maclaurin spheroids, as can be 
found in Chandrasckhar (1967a) (C67a) equation 
(3), where is a constant containing a Newtonian 
and post-Newtonian contribution, cf. (28) of that 
paper. 

Next we note, that for a given equation of state, 
an axially symmetric, stationary and uniformly ro- 
tating fluid is described by two parameters. For 
our purposes, we can take them to be a^/ai, which 
we prescribe using (4) , and the value for ai , which 
we leave undetermined. 

One has two additional degrees of freedom, 
which amount to the mapping between a Newto- 
nian and post-Newtonian solution and is a matter 
of convention (cf. Bardeen 1971). For example, 
one can write the coordinate volume of the star to 
be 



V = Vo + Vi5 + .. 



(31) 



where 5 is some relativistic parameter, and then 
choose to have the PN contribution vanish, Vi = 0. 



This is the choice that was made in CE78 and 
C67b and also in Chandrasekhar's original paper 
on the PN Maclaurin spheroids C67a. We have 
followed this covention in the current paper, mak- 
ing it easy to compare our results to those of C67a. 
The second degree of freedom one has was left un- 
specified in much of C67a, though Table I lists 
values with the choice Sf^ = = 0.^ 
If we introduce the new coordinate 



(32) 



and make use of (13), (15) and (17), then the 
bounding surface (cf. (47) in CE78) is given by 



= — + 4 

5 
3 

Ualxj 
9 ai 



2a\x\ 



+ 5: 



5 2 



TJJ 



3^^ 

4-1 

ai 



■ 4^7 



2-^3 



(33) 



Using the equation for the surface Tu'^/af = 1 — 
3:3/03, which holds at the Newtonian level and can 
thus be inserted into the PN term above, one sees 
that the term with xf vanishes and one finds that 
the equation is identical to (42) of C67a if 



13 



:2'^3 



13a^ 



and 



C cM , l^ oM 
bi= + —b2 - 



13of 



(34) 
(35) 



hold. As mentioned in that paper, 5*3 = may 
be chosen without loss of generality^ which then 
leads to a unique relationship between ^3 and 6*2 ' , 
which is shown to be correct in Appendix B3. The 
constant Sf^ can be chosen arbitrarily just as with 
(which depends on wi). 

If one considers the limit 02 — >■ ai and simul- 
taneously requires that the star rotate uniformly. 



® Where necessary, we distinguish the constants of C67a from 

those used here by adding the superscript 'M'. 
'^Note that (34) and (35) together with (15) are equivalent 
to the three equations (123) of C67b as can be seen either 
or identifying a of that equation with 



by taking 5f 



5f 



and /3 with 



M °-3 qM 



5 



then (28) provides the unique value for wi, 



wi w 0.01646051799, (36) 

which is equivalent to making a choice for Sf^ dif- 
ferent from the one made in C67a, but no more 
and no less physically meaningful. 

The most significant result of the analysis of 
the axisymmetric limit is that (28) shows us that 
the rigidly rotating limit (ri = 0) and the original 
choice of velocity field in CE78 {wi = W2 = 0) 
are incompatible. While it is possible with that 
velocity field to find the post-Newtonian A/[aclau- 
rin solution at the bifurcation point, this solution 
is not continuously connected to any other solu- 
tion. When considering the question of the ex- 
istence or non-existence of non-axially symmetric 
but stationary solutions, it seems important to re- 
tain the possibility of studying a neighbourhood 
of the axially symmetric and uniformly rotating 
limit, especially since such solutions are known to 
exist^. This possibility was excluded by the ap- 
proach taken in CE78. 

In a follow-up paper, we intend to tackle 
the problem with a more general approach that 
lends itself better to proceeding to higher post- 
Newtonian orders, is not as restrictive in the so- 
lutions it permits and allows one to show that the 
singularity discussed in CE78 is an artefact of the 
specific method chosen and not an inherent prop- 
erty of the post- Newtonian Dedekind solutions (cf. 
Giirlebeck & Petroflf 2010). 

We gratefully acknowledge helpful discussions 
with M. Ansorg, ,J. Bicak, ,J. Friedman and R. 
Meinel. The first author was financially supported 
by the grants GAUK 116-10/258025 and GACR 
205/09/H033 and the second by the Deutsche 
Forschungsgemeinschaft as part of the project 
"Gravitational Wave Astronomy" (SFB/TR7- 
Bl). 



'As far as we know, there exists no formal proof demonstrat- 
ing the existence of such solutions. Steps in that direction 
were taken by Hcilig (1995) and the existence has been 
demonstrated by many groups that arc able to solve Ein- 
stein's equations numerically to extremely high accuracy, 
see e.g. Ansorg et al. (2003) 
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A. A Detailed Discussion of Chemdrasekheir and Elbert's Work 



Wc mentioned in footnote 3 that we have been unable to reproduee the values from Table 1' in CE78 
that result from solving (58') nor have we succeeded in finding the source of the discrepancy. It is important 
to rule out an error in our understanding of that paper or an error in our own solutions to the equations 
presented there, and we therefore provide a detailed discussion here (in this section we use the velocity field 

in that paper, i.e. wi = W2 = 0). 

The calculations we performed were done with the aid of computer algebra. As a test, we did all the 
calculations using both Maple and Mathematica. To be absolutely certain that we solved the equations 
correctly, we wrote down the line element and energy-momentum tensor as given in Chandrasekhar (1965), 
had Mathematica (TTC package) determine Einstein's equations to first post-Newtonian order and then 
verificid that they are indeed fulfilled. When the values from Table 1' of CE78 are inserted, then one finds 
that the condition that pressure vanish on the surface is violated at a level three orders of magnitude higher 
than with the values from our Table 1. We also verified that the violation vanishes in our case as more 
significant figures are added. 

The solutions we found for gi, q2 and agree with those given in Table 1' of CE78. This provides strong 
evidence suggesting that our numerical evaluation of a^/ai for a given aijax and of the index symbols is 
correct. Moreover, the dependence of q, ri, r2, t\ and ti on Si as given in equations (37') and (53')-(56') 
can be seen to hold both in Table 1' and Table 1. This indicates strongly that the typo in equation (53') 
of CE78 mentioned in footnote 5 is truly only that and that the quantities in the integrability condition of 
(11') are treated correctly in both papers, leaving only W and $ to be verified. 

The system of linear equations providing the values for 5,, i.e. (58'), can of course be written as follows: 



For a given value of a^jai^ the matrix (J^lij) depends on the m's from C67b and via their Si dependence, 
indirectly on g, ri, r2, t\ and ti- The vector (iVj) depends on the a's and again on the (non-S'j dependent 
part of) g, ri, r2, t\ and t^. We return to a discussion of this equation after mentioning a few incongruities 



In (44') a factor 1/{ttGp) is missing in 5U because the equation is copied directly from (74) of C67b, 
whereas the relationship between p/ p and 5U is not the same in (39') of CE78 and (75) of C67b. This 
mistake is corrected in (45') and (46') however. In (39') there is also a factor l/(7rGp)^ missing in the term 



2$ + 2v^{7 + M - j as can be seen by checking dimensions^ and comparing to (11) in Chandrasekhar & 



Elbert (1974). Finally, we note that (Al) from above only ensures that the pressure is constant on the 
surface of the PN-ellipsoid as discussed in C67b, cf. (75) in loc. cit., but not that it vanishes. The constant 
that would have to be determined to ensure vanishing pressure was not written in (39') or (40')^'^ and the 
constant that is a part of SU in (44') was dropped when proceeding to (45'). Since the determination of this 
constant plays no role in the paper however, we need not discuss it further and have not done so in our own 



We find that the determinant of (My ) vanishes at a^jai = 0.33700003168 . . . just as in CE78, where it 
is given to four significant figures. This provides evidence suggesting that the matrices agree (and thus the 
6U) and that the vectors {Ni) disagree. If we multiply 6U by a factor tt, as suggested in the last paragraph. 



'We advise the reader that, as mentioned after (14) in Chandrasekhar & Elbert (1974), the units in which Qi and Q2 are 
measured change as of this point by a factor sJ-KGp. 

'The constant contained in W is completely determined by (44') and is thus not available as a variable to ensure that the 
pressure vanish on the surfax;e. 




(Al) 



in GE78. 




paper. 
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then the determinant becomes zero for a-ijax = 0.30874 Nonetheless, we tested that neither an arbitrary 

factor in front of this term, nor one in front of the term 2$ + 2v^U + | ^ can explain the results in CE78. 

A natural explanation for a disagreement between the vectors (Ni) in our case and in CE78 would be 
that one of the a's contains a mistake. We checked to see that an arbitrary change in a single a cannot 
account for the differences in the results however. Since an explicit expression for these a's is not provided 
in CE78, we cannot test directly to see whether or not each agrees. However, in the implicit expressions 



2 



from (39') and (40'), only the contributions from 2$ + 2v^U + \ j are not written out. These can easily 
be compared to those written out explicitly for the a's of C67b, where the appropriate modifications for the 
different Newtonian velocity have to be taken into account, and show perfect agreement with our expressions. 
In particular, the relationship to the a's of C67b for 02 = oi, which is discussed in Appendix A.l provides 
additional evidence for the correctness of our expressions. We also generated the a's with computer algebra 
by typing out the expressions for (11'), solving the integrability condition and integrating it and showed that 
these agree with the expressions provided below. 
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af = - AaiBi2{Ai + A2) - 2IQ1Q2 - {21 + aa^As)^! 



a 



(A2) 



(A3) 
(A4) 



+ {afAu - ^Bn)i2alQl - 2a\A^ - ZalA^) - \B^2{^alQ\ - 2a^^2 - Ba^^s) + lalA^Bxz 

al^ = — ^Aj - 4alBi2{Ai + A2) - 2IQ1Q2 - {21 + 3ai^3)^2 

+ {alA22 - lB22){2alQl - 2a\A2 - 2,alA^) - iSi2(2o?Qi - 2a\A^ - ZalA^) + fa^AsBsa 
af = -a\Ai - {21 + 2.alAi)Ai - \B23{2alQ\ - 2a^^2 - Ba^^s) 

- \Bis{2a\Ql - 2a\Ai - ialA^) + |a2^3(B33 - 20^^33) 

all = -?2^3 - 2Q? (ai + ^^2) + {2a\Ql - 2af^i - 3a^A3)(-o?An2 + 5B112) 

+ {2alQl - 2alA2 - 3a§^3)(-aiAi22 + 5^122) - iai^3Si23 + 2Q1Q2 - ^) A2 (A5) 

- \QlQ2 - 2aiQi02(3^22 + A12) + AQl (^A^ - ^^A^ - Qi{qi + ^q2) 

all = 4^1 - 2(3?^3 + 2Q? (1 - 4 V3 - 2alQiQ2{Ai3 + A23) + {2a\Ql - 2a\A^ - 2,alAi){\B^2i) 
02 \ 

+ {2alQ\ - 2alA2 - 3a^A3)(-aiA223 + 5-B223) - f a^A3(-2a^^233 + ^233) 

(A6) 

all = ^Al - 2QlAs + 2Q1Q2 f 1 - ^ ^3 - ^alQMA.a + A23) - |a^A3(-2a^Ai33 + B133) 

a-i \ a^y (A7) 

+ {2alQl - 2ajAi - 3a2A3)(-a?An3 + i^na) + (2a^0? - 202^2 - 3alA3)(iBi23) 
all = k^^l - '^Ql^^ + (2a?<32 - 2a?Ai - 3a^A3)(-a2>lni + ^B^ii) 

+ {2a\Q\ - 2alA2 - ZalAs){\B^i2) - ^alA^B^^^ + Q^Q2 [l - ^) A^ ^^^^ 

- \QiQl - a2<3iQ2(3Aii + A12) - \Q2q1 



For the terms in (11'), we checked our expressions by ensuring that V^J7 = —inGp, (8') and the Newtonian equations hold. 
Furthermore, we tested the u's by first ensuring that the moments S,, Dijk fulfil the appropriate Poisson equation and that 
the SU^"^^ of (69) from C67b agree with (70) and (71) from the same paper. 
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all = 14^3 - 2Q?A2 + i2alQl - 2af Ai - 3a^A3)(iSi22) 

CI') 



+ {2alQl - 2alA2 - Sa^ A3) (-0^^222 + ^-8222) - 103^3^223 + (^1 - ^) ^2 



(A9) 



- iQ?<32 - QiQ2af (3^22 + A12) - iQig2 
^alA3{-2alA333 + ^^333) 



4'' 

^3 = 1^3 + 3(2a?Q2 - 2a? - 802^3)5133 + i(2aiQ2 _ 2alA2 - 80^^3)5233 

5„2 . / « 2 . . 1 „ X y^^^) 



Let us summarize the arguments from above. We have checked all the equations in Part I of CE78 and 
find the analytic expressions to be free of error, except for the few minor points mentioned above. We 
have good reason to believe that both in that paper and here, Einstein's PN-equations are solved correctly 
including the PN-Bianchi identity. We obtain different numerical values for Si which we suspect is related to 
a problem with the numerical evaluation of the a's in CE78, though we cannot be certain that our matrices 
{Mij) agree simply because their determinants vanish at the same point. The various tests of our a's and 
the fact that we find the post-Newtonian Maclaurin spheroids in the axisymmetric case convince us that our 
values are correct. 

A.l. The Solution at the Bifurcation Point 

At the point 02 = Oi, i.e. at the bifurcation point along the Maclaurin sequence, the following relations 
can be used to simplify the expressions for the a's, where Q refers to the angular velocity of the uniformly 
rotating Newtonian solution and has the same meaning as in C67b: 

alAs = ajAii = al{Ai - Bn) = I - 2alAi = al {Ai - ^Ql) , = 2Bii = Qj. (All) 

Note that at this point, the a's of C67b and C67a agree and we find 

al^ = - SafBnAi + 2/0? - {21 + 802^3)^1 + (a?ylii - Sii)(2a?0? - 2a?Ai - 802^3) + |a3^3Bi3 

Oi 



-15o?A? - falQi + UafAiQl + ^ajAsB 
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= af = a^s = ^67 



(A12) 



al^ = -alAl - {21 + 8o^A3)A3 - Bi3(2af - 2o?Ai - SajAs) + §0^^3(533 - 201^33) 

2/r,^2 ir, A \ A , „4^2 



= ai {2Qi - 10^1)^3 + atQiAis (A13) 



:af 



a* 



all = -i^l - 4Q?^i + aj{lQl - 5Ai){-2alAni + Bin) - |a^^3Sii3 + ^Qt + 8afQf An - ^iQi 



= a 



(A14) 



alt = - 20?^3 + 4a?g?Ai3 + a?(|g? - 5Ai)(Bn3 - o?Aii3) - fo^ A3 (-20^^133 + ^133) , , ^ 

Oi (Aloj 

^67 ^78 „67 
= "23 = "31 = "31 

all = all = Wil + 5^191 = afl + \Qm (A16) 



41 = 1^3 + WiilQl - 5Ai)Bi33 - fa2A3(-2a2A333 + 1^333) 



= a 



(A17) 
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Table 1 

The numerical values we find for the quantities listed in Table 1 of Chandrasekhar & Elbert (1978). 
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T2 






1.00 


2, 


.8272 


-2, 


.8272 





.0000 


-0.3050 


-0.3290 


-0.0144 


0.0574 


0.2398 


-2.7984 


0.0288 


-0.0288 


0.0000 


0.0000 


0.99 


2, 


.8173 


-2, 


.8370 


0, 


.0211 


-0.2944 


-0.3323 


-0.0132 


0.0542 


0.2481 


-2.7984 


0.0196 


-0.0378 


-0.0454 


-0.0445 


0.98 


2, 


.8073 


-2, 


.8470 


0, 


.0424 


-0.2838 


-0.3355 


-0.0120 


0.0508 


0.2565 


-2.7984 


0.0101 


-0.0466 


-0.0923 


-0.0887 


0.97 


2, 


.7972 


-2, 


.8570 





.0639 


-0.2733 


-0.3387 


-0.0107 


0.0474 


0.2648 


-2.7984 


0.0002 


-0.0553 


-0.1407 


-0.1324 


0.96 


2, 


.7869 


-2, 


,8671 





,0857 


-0.2628 


-0..3416 


-0.0093 


0.0438 


0.2732 


-2.7984 


-0.0101 


-0.06,38 


-0.1906 


-0.1757 


0.95 


2, 


.7766 


-2, 


,8774 





,1077 


-0.2524 


-0.3445 


-0.0078 


0.0401 


0.2815 


-2.7985 


-0.0207 


-0.0722 


-0.2422 


-0.2186 


0.90 


2, 


.7231 


-2, 


,9297 


0, 


,2213 


-0.2010 


-0.3573 


0.0010 


0.0200 


0.3223 


-2.7989 


-0.0804 


-0.1121 


-0.5287 


-0.4282 


0.85 


2, 


.6665 


-2, 


,9843 


0, 


,3414 


-0.1509 


-0.3675 


0.0121 


-0.0030 


0.3610 


-2.7999 


-0.1528 


-0.1486 


-0.8741 


-0.6316 


0.80 


2, 


.6067 


-3, 


,0412 


0, 


,4691 


-0.1019 


-0.3754 


0.0258 


-0.0288 


0.3949 


-2.8020 


-0.2413 


-0.1818 


-1.2996 


-0.8317 


0.75 


2, 


.5439 


-3, 


,1006 


0, 


,6052 


-0.0538 


-0..3813 


0.0426 


-0.0572 


0.4198 


-2.8064 


-0.3506 


-0.2120 


-1.8363 


-1.0329 


0.70 


2, 


.4781 


-3, 


,1627 


0, 


,7511 


-0.0064 


-0..3856 


0.0636 


-0.0871 


0.4280 


-2.8152 


-0.4872 


-0.2398 


-2.5324 


-1.2409 


0.66 


2, 


.42,36 


-3, 


,2146 


0, 


,8758 


0.0313 


-0..3882 


0.0845 


-0.1106 


0.4131 


-2.8280 


-0.6222 


-0.2610 


-3.2547 


-1.4177 


0.65 


2, 


.4098 


-3, 


,2279 


0, 


,9082 


0.0408 


-0..3888 


0.0905 


-0.1161 


0.4049 


-2.8324 


-0.6603 


-0.2662 


-3.4656 


-1.4642 


0.64 


2, 


.3958 


-3, 


,2414 





,9412 


0.0503 


-0..3895 


0.0968 


-0.1215 


0.3944 


-2.8374 


-0.7005 


-0.2715 


-3.6910 


-1.5118 


0.63 


2, 


.3818 


-3, 


,2550 


0, 


,9747 


0.0598 


-0.3901 


0.1037 


-0.1265 


0.3814 


-2.8432 


-0.7427 


-0.2768 


-3.9323 


-1.5607 


0.62 


2, 


.3677 


-3, 


.2687 


1, 


.0087 


0.0693 


-0.3907 


0.1109 


-0.1311 


0.3653 


-2.8497 


-0.7873 


-0.2821 


-4.1913 


-1.6111 


0.61 


2, 


.3536 


-3, 


.2826 


1, 


.0434 


0.0789 


-0.3914 


0.1187 


-0.1352 


0.3459 


-2.8573 


-0.8343 


-0.2876 


-4.4697 


-1.6632 


0.60 


2, 


.3394 


-3, 


.2967 


1, 


.0786 


0.0886 


-0.3921 


0.1271 


-0.1388 


0.3225 


-2.8659 


-0.8840 


-0.2931 


-4.7697 


-1.7171 


0.59 


2, 


.3252 


-3, 


.3109 


1, 


.1145 


0.0983 


-0.3929 


0.1362 


-0.1416 


0.2947 


-2.8758 


-0.9367 


-0.2989 


-5.0939 


-1.7732 


0.58 


2, 


.3109 


-3, 


.3254 


1, 


.1510 


0.1082 


-0.3938 


0.1460 


-0.1436 


0.2616 


-2.8872 


-0.9924 


-0.3048 


-5.4449 


-1.8317 


0.57 


2, 


.2967 


-3, 


.3400 


1, 


.1881 


0.1182 


-0.3948 


0.1567 


-0.1444 


0.2224 


-2.9004 


-1.0516 


-0.3110 


-5.8262 


-1.8929 


0.56 


2, 


.2824 


-3, 


.3548 


1, 


.2260 


0.1284 


-0.3960 


0.1685 


-0.1439 


0.1763 


-2.9155 


-1.1146 


-0.3175 


-6.2416 


-1.9574 


0.55 


2, 


.2681 


-3 


,3699 


1, 


.2646 


0.1387 


-0.3973 


0.1814 


-0.1418 


0.1219 


-2.9330 


-1.1818 


-0.3244 


-6.6957 


-2.0254 


0.50 


2, 


.1975 


-3, 


,4487 


1 


,4690 


0.19-56 


-0.4087 


0.2720 


-0.0912 


-0.3,361 


-3.0751 


-1.5988 


-0.,3691 


-9.7775 


-2.4444 


0.45 


2, 


.1301 


-3, 


,5353 


1 


,6959 


0.2722 


-0.4377 


0.4537 


0.1171 


-1.4497 


-3.4139 


-2.2489 


-0.4521 


-15.4156 


-3.1217 


0.40 


2, 


.0694 


-3, 


,6327 


1, 


,9502 


0.4253 


-0.5330 


0.9625 


0.9405 


-4.8952 


-4.4794 


-3.5956 


-0.6789 


-29.1295 


-4.6607 


0.35 


2, 


.0208 


-3, 


,7456 


2, 


,2386 


1.6155 


-1.4902 


5.7347 


10.0089 


-37.9506 


-15.1294 


-13.4666 


-2.7757 


-137.5070 


-16.8446 


0.34 


2, 


.0132 


-3, 


,7707 


2, 


,3012 


6.5850 


-5.5917 


26.0552 


49.3174 


-178.2659 


-60.8513 


-54.2058 


-11.6728 


-582.0542 


-67.2855 


0.33 


2, 


.0065 


-3, 


,7968 


2, 


,3655 


-2.6453 


2.0383 


-11.7342 


-23.8765 


82.4680 


24.2357 


21.4267 


4.8654 


241.2314 


26.2701 


0.32 


2, 


.0007 


-3, 


,8240 


2, 


,4317 


-1.0183 


0.6975 


-5.0896 


-11.0416 


36.4961 


9.3003 


8.0964 


1.9536 


94.8022 


9.7077 


0.30 


1, 


.9925 


-3, 


,8819 


2, 


,5703 


-0.4078 


0.1986 


-2.6145 


-6.2981 


19.0909 


3.7752 


3.1511 


0.8615 


37.3,393 


3.3605 


0.28 


1, 


.9892 


-3, 


,9456 


2, 


,7179 


-0.2340 


0.0578 


-1.9180 


-4.9749 


13.7746 


2.2579 


1.8969 


0.5449 


17.7005 


1.3877 


0.25 


1, 


.9960 


-4, 


.0540 


2, 


.9584 


-0.1417 


-0.0205 


-1.5433 


-4.1994 


9.8393 


1.4994 


1.8182 


0.3549 


0.1194 


0.0075 



B. Explicit expressions for Si, S3 and ri 

At the point 02 = ai, the u's from C67b and C67a are related by 



(2)M 



_9_ 
13 



,(3) 



4u, 



(4) 



17a| (5) 
3ai ""'^ 



(Bl) 



which follows from (119) of C67b. 

Using these relations, those between the a's and (11), (13) and (17), one finds that the third of equations 

(58) of CE78 becomes 



= atQll^alalQll + atQll 

a 4 67 2 2 67 , 4 67 , ^l*3l , ^^O 2 2 a o ^'^ / 4 (2)M 2 2 (2)M , 4 (2)M\ ,-r,ri\ 

= atali - aia^a% + a^a% + + —afa^AaSa - — {atu\i' - afa^uy + aluy j S3 (B2) 

We thus have the solution 

9 afall - alalall + a^atl + aiQin/2 



S3 = 



(B3) 



which agrees with (99) of C67a if we take (34) of this paper into account. In order to provide concise explicit 
formulae, we again make use of the eccentricity 



the quantity 



and recall that Qi is 



C := 1046*^ - 444e'' + 6306^ - 245 



(B4) 



8e2(l - e2) 

3 + 8e2 - 8e4 



We now provide explicit expressions for ^i, 6*3 and ri. Note that the expressions for 5*1 and 6*3 can be 
obtained either as limiting values or by placing oneself directly on the point 02 = ai . On the other hand, ri 
can only be obtained by a limiting process. The formulae read 



Si 



2e2 



-1 



26eC 



(2864e* - 10128e^ + 14712e^ - 81206^ + 1365)Q? 



e 4e 



(224e'' - 840e^ + 1170e^ - 455)ri 



a 366^ 



65C 



(272e^ - 2446^ + 35)Q? 3e 



n 



8e2 

-(24e^ - 12e^ - 1) 



Qi 
7 
4 



Wi. 



8e2(2e2 + l) 

In deriving these expressions, we have made use of the identities 



3^? - 3^1 - 4a?^iAii + ?>a\Aii - 2a\A\A 



= 0, 



0. 



(B5) 

(B6) 
(B7) 

(B8) 
(B9) 
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